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Abstract 

We calculate the time evolution of the entanglement entropy in a 1+1 CFT with a holographic 
dual when there is a localized left-moving packet of energy density. We find the gravity result agrees 

(N 

with a held theory result derived from the transformation properties of Renyi entropy. We are able 

Oh 

to reproduce behavior which qualitatively agrees with CFT results of entanglement entropy of a 
system subjected to a local quench. In doing so we construct a hnite diffeomorphism which tales 
three-dimensional anti-de Sitter space in the Poincare patch to a general solution, generalizing the 

-t— > 

diffeomorphism that takes the Poincare patch a BTZ black hole. We briefly discuss the calculation 

<D 

& of correlation functions in these backgrounds and give results at large operator dimension. 
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I. INTRODUCTION 

The gauge-gravity correspondence [1] tells us that certain quantum field theories have 
dual descriptions in terms of a quantum theory of gravity. Because this is a weak-strong 
duality, we can use semiclassical gravity to learn about strongly interacting field theories, and 
we can connect black hole thermodynamics with thermodynamics of the dual field theory. In 
the limit where we have control over the gravitational side, we can use the correspondence to 
understand nontrivial behavior of systems at strong coupling. One thing in particular we can 
use this to study is the behavior of the system in nontrivial and possibly time- dependent 
excited states. By studying these states we can gain insight about the dynamics of the 
strongly interacting theory. 

One nontrivial way of probing the state is by studying the entanglement entropy of a 
spacial region A with it's complement B = A. This is defined by the von Neumann entropy 
of the reduced density matrix Sa = — tr(p^ lnp^), where pa = tr^p. The conjectured 
holographic prescription for calculating entanglement entropy for a static configuration [2] 
involves finding minimal surfaces E in the bulk geometry whose boundary coincides with 
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the boundary of A, OA = <9£. The entanglement entropy is then 1 

S A = ^-mm[Vo\(E)} (1.1) 

This is only a conjectured formula, though it has passed many nontrivial checks (see e.g. 0,11] 
for a review). When the system in question is no longer static, this must be generalized to 
the more nontrivial covariant holographic entanglement formula [5]. However, it was argued 
that in three-dimensional gravity without extra matter fields, the fact that any solution is 
in fact Ad S3 (up to quotients) greatly simplifies calculating the entanglement entropy. We 
will end up using this fact to greatly simplify our calculations. 

The goal of this paper is to calculate the entanglement entropy of a conformal field theory 
as a pulse of energy passes through it by using holographic techniques. We demonstrate the 
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FIG. 1: A representation of the system we wish to construct. The red shaded diagonal area is the 
region of nonzero left-moving energy density with a (lightcone coordinate) width A, and the green 
vertical region is the time-swept path of the spacial regions we wish to calculate entanglement 
entropy for. 

construction in figure [I] This problem is significantly more tractable in AdSz/CFT2 by 
taking advantage of the Weyl anomaly in the stress tensor transformation. It is also com- 
parable to the time dependent entanglement entropy in a CFT subjected to a local quench 

1 We are working in Einstein frame with the action normalized as S = 167r 1 Gw / y/—gR + . . . and ignoring 
higher-derivative corrections. 
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[61 El- In section [fi] we will review the enhanced conformal symmetry in asymptotically AdS% 
gravity. We will construct a finite diffeomorphism that takes us from the Poincare patch to 



a general solution with arbitrary T±±. In section III we use this diffeomorphism to construct 
a solution corresponding to a left-moving packet of energy with compact support. Using the 
explicit diffeomorphism we are able to calculate the time dependent entanglement entropy of 



a region as the wavepacket passes through it. In section [TV] we take the narrow width limit 
of the wavepacket, reducing it to a three-dimensional Aichelburg-Sexl shockwave, and argue 
that this result is related to the time evolution of entanglement entropy after a local quench. 
In section [V] we find that the transformation law derived holographically can also be found 
by field theory considerations, giving a nontrivial verification of the covariant holographic 



entanglement entropy formulation [5J. In section VI we discuss implications and various 
generalizations of our results. 

II. 1+1 CONFORMAL SYMMETRY AND AdS 3 

Consider general relativity in 2+1 dimensions with a negative cosmological constant in 
a spacetime that's asymptotically AdS%. Thanks to the fact that there are zero degrees 
of freedom in 2+1 GR, the most general solution to Einstein's equations is (up to gauge 
transformations) [S] 2 

ds 2 = f ^L+dxl + L_dxl ~(J^ + dx + dx - + ^r) (2- 1 ) 

where L± = L±(?/±) and the holographic stress tensor [5] is 

T *± = »k L±,T * = - (2 - 2) 

Of course this manifold is still (locally) exact AdS^, we have just chosen a nontrivial slicing. 
We can in fact shift L± with a linearized diffeomorphism, 

/ ~2fl -Aj^ ill \ / fl _|_ fl \ 

which shifts L± as 

rill 

AL ± = f ± L' ± + 2L ± f' ± - J -^-. (2.4) 

2 The paper of Banados works primarily in Euclidean signature but it is simple to extend the results to 
Lorentzian signature. 
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Integrating this diffeomorphism to nonlinear order in a closed form is not trivial, though 
it has been done for the simple case of a constant L±, which is just the change of slicing 
to go from Ad S3 to the BTZ black hole [ID] . We find that it is also possible to construct 
the nonlinear diffeomorphism when the starting coordinate slicing is that of Poincare AdS?,. 
Starting with the metric 

ds 2 = f(-2dy + dy. + du 2 )/u 2 , (2.5) 
acting with the coordinate transformation (y + ,y_,y) — » (x + ,X-,z) given from 



2z 2 r 2 r; (4/;/:) 3 /2 

— — ~ 11 = z ' — 



V± = f±(x±) + - 7 , 7 , ± T , u = z + (2.6) 



which linearizes to (2.3) with the original L±(y±) = 0, and asymptotes to 

V± = /±(x±) + 0(z 2 ), u = z^ffJL + 0(z 3 ), (2.7) 



the metric is precisely (2.1) with L±{x±) proportional to the Schwarzian derivative of f. 



2^' ^ 4/ : 



We see from (2.2) the stress tensor transformation satisfies the general transformation law 
with T a b(y±) = 0. When the theory is on a cylinder, we have to be very careful about how 
the map changes the boundary conditions, so we will consider the theory on Ai 1,1 . Thanks 
to our closed form diffeomorphism we know solving any problem such as wave scattering 



or geodesic motion on (2.1) reduces to solving that same problem on Poincare AdS and 



then following the coordinate map, often paying careful attention to what we do to the 



boundary. We can use (2.6) to map solutions to, for instance, the Klein-gordon equation 
in y coordinates to nontrivial x coordinates. However, as we are working in Lorentzian 
signature, there will be nontrivial subtleties in terms of boundary conditions on the horizon 
if we are interested in calculating different-time correlators. Thanks to the closed-form finite 



diffeomorphism these can be sorted out by inspection once we solve (2.8). 

It needs to be emphasized that there can be other states in the CFT with the same stress 
tensor expectation values, which may differ by the behavior of other operators. In the bulk 
this would correspond to turning on other supergravity fields. These other solutions will 
not be related to planar AdS by a simple diffeomorphism, and so we can not as easily study 
entanglement in these backgrounds. However, when the scale of the stress tensor expectation 
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value is much greater than the one-point functions for other operators, this background will 
be a good approximation. 3 



III. LEFT-MOVING "WAVEPACKETS" 



Let's start by constructing a simple left-moving wavepacket with a constant nonzero T ++ 
for < x + < A, and L_ = everywhere. We can think of this as an AdS^ version of an 
AdS Aichelburg-Sexl shockwave [H] , where the stress tensor is spread over a finite region of 
x + and not just delta-function supported. The simplest model is one where 



| x + < 

r 2 | < x + < A , L. 

1 A < x+ 



0. 



(3.1) 



This is only piecewise smooth, but we know that we could easily replace L + with a function 
that is smooth and still has compact support. We will find that this discontinuity in L + 
means that the entanglement entropy has a discontinuous second derivative in x + . We now 
simply need to solve ( |2.8 ) with the source given by (3.1 ). First, let's point out the diffs that 
give L + = are /+ = (a + bc + )/(c + dc + ) and similarly for the right-moving sector 4 . Because 
we are not exciting the right-moving sector, we can happily pick /_ = x_. For x + < we 
will simply use /+ = x + and a general SL(2, R) element for x + > A. What about a positive 
constant L + = r 2 ? A general solution is 



/+ = a + b tanh [tx + + c] 



(3.2) 



We need T ++ to be piecewise continuous, so matching / + , f' + , /" at x + = and x + = A 
we find 

fx + x + < 
tanh[rx+]/r | < x + < A (3.3) 

x+ > A 



1 l + (x + -A)rcoth[rA] 



t (x+-A)T+coth[TA] 

We demonstrate this map in figure [2j Note that the map does not cover all of y + , only 
y + < coth[rA]/r. This is similar to what happens if we act with a diffeomorphism that 



3 We thank Matthew Headrick for emphasizing this to us. 

4 This is just the global SX(2,R) x SX(2,R) symmetry. 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

FIG. 2: The map /+ for A = 1/4, r = 2. Note that the map is the identity for x+ < 0. On the left 
we also show the region with nonzero T ++ and on the right we shade the region of x not covered 
by the y coordinates. 

takes us to an extreme left-moving BTZ black hole, where /+ = exp(-y/L + x + ). We know 
that geodesies in y coordinates connecting two points on opposite sides of this boundary 
must map to a geodesic which begins on the x boundary and crosses a horizon in the 
interior [12J 5 . However, since /+ is invertible a geodesic connecting two spacelike-separated 
boundary points in x coordinates will map to a spacelike geodesic connecting two boundary 
points with both y + < coth[rA]/r. We can demonstrate this by using / + , its inverse, and 
boosts to map geodesies connecting points at equal y± time, u(y) = \/(y — Vi){y2 ~ y)i to 
geodesies connecting two equal x± time points in x coordinates. The explicit functional form 
of the geodesies is messy and not terribly enlightening, so we provide plots to help visualize 
in figure [3j 

Now we wish to measure the entanglement entropy for a region A between two points 
in x±,z coordinates. This is related to the (regulated) AdS geodesic distance between the 
two points by simply mapping the problem to y coordinates. This is quite simple, as in y± 

5 It may seem as though acting with a global SL(2, K) rotation might alleviate this, but since an element 
of SL(2,R) maps MlJ{oo} H> MlJ{oo}, we know even with this rotation the range of this function will 
never cover all K. 
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FIG. 3: Geodesies connecting two points at equal x time, with r = 2, A = 1, d = 3. On the 
left we have the projection to the (x, t) plane, with dashed red lines denoting where the region of 
nonzero L + is. On the right we have the projection onto the (x, z) plane, with z shifted by time 
for clarity. Note that while the curves look timelike, this is due to the frame-dragging from L + 
and they are indeed spacelike. Note that when A > d the curves are very similar, and when they 
are entirely within the pulse they are simply extermal BTZ geodesies. 



coordinates the answer (for spacelike seperated points) is 



Sent. (?/l,?/2) = ^Og 



4(2/2 



2/1- 



) 2 (2/2- 



2/1- 



2 2 

fc 2 fc l 



c = U/2G 



N ■ 



(3.4) 



where q is the location of the radial cutoff assumed to be much smaller than any other 
bulk length scale. It is important to keep track of these, as our answer in terms of x 
coordinates is meant to be asked with the same regulator on each point in z, not u. Note 
that if we calculated an object which is UV finite such as a mutual information, these terms 



will cancel but it is possible that subleading terms in (2.6) may be relevant, which would 



indicate some regular ization-scheme dependence of the result. 



Let's first work out how (3.4) transforms under a general left-moving transformation / 4 



with /_(x_) = x_. The result is 



S en t.\ x l: x 2) 



12 



log 



4(/+(x 2 +) - /+(x 1+ )) 2 (x 2 , - xx-l 
f' + {x 2+ )f' + {x 1+ )5l8l 



(3.5) 
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where 5i is the radial cutoff in z% and we have used the fact that we require 



5 <C (scales in / + ). 



(3.6) 



This seems like it should generalize when /_ is nontrivial as well to 

"4[/+(x 2+ ) - / + (s 1+ )] a [/-(*2-) - f-(xi-)? 



S e nt.\ x li x %) 
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log 



(3.7) 



f^x 2+ )f + (x 1+ )fL(x^)f_(x^)6 2 5 2 

which would be true as long as we can verify that the connected geodesic between the two 
points is always shorter than a piecewise smooth geodesic which falls down into the horizon 
and runs along it. However we will find from CFT considerations that this is precisely the 
behavior we find. Recall that when g e SX(2,M), 

(g( x ) -g(y)) 2 



[X 



yf 



(3.* 



g'(x)g'(y) 

and so despite swirling the coordinates around, the regulated geodesic distance between two 
boundary points is invariant. This means that our entanglement entropy will go back to 
the vacuum value when the pulse has completely passed through the region, despite the 



coordinate transformation / + no longer being simply the identity. Indeed, plugging (3.3) 



into (3.5) and evaluating at equal times, 



t\ = t 2 , Xi = x 2 + V2d, i.e. X\± = x 2 ± ± d 



we find 



S A =°- log [2D 2 /5 2 ] , AS A = S A - S v 
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\og[D 2 /d 2 



(3.9) 



(3.10) 



D 



d 2 

d cosh [rx + ](d + tanh[rx + ]/r — x + ) 
dwah.\rd]/r 
c?cosh[rA] 

x (d - A + ( 1 +( d ~ x +)( x +- A ) r2 ) tanh[rA] 

dcosh[r(d + A — x + )] 
x ( x+ - A + tanh[r(d + A ~ x+)1 
d 2 



x + <0 (I) 
< x + < A, x + < d (II) 
0<x + <A,0<x + -d<A (Ilia) 
A < d < (II lb) 

A < x+, A < x + -d < (IV) 



A + d < x+ (V) 



(3.11) 
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FIG. 4: The six regions for the piecewise smooth function (3.11), labelled I - V. The region of 



interest is always between x + — d and x + . Note that Ilia only occurs for d < A and Illb for A < d. 



To clarify this we include a figure showing what regions we are considering in figure |4j Note 
that .D 4 is symmetric under a reflection about x + 
gives us the peak entanglement entropy 



^4^, where it has a maximum. This 



max(AS') 



log 



sinh [rd] 



log 



A(d- A)r cosh[r A] + (4+(rf- A) 2 t 2 ) sinhfrA] 
4-rd 



d < A 
A < d 



(3.12) 



We also provide plots of AS for various cases in figure [5] We again emphasize that this is 



the entanglement entropy for the state created by acting with the diffeomorphism (3.3) on 



the vacuum, and it is possible that there are other states with the same stress tensor profiles 
but with for instance nonvanishing one-point functions for other operators. However for 
states where the stress tensor scale is much greater than any other scales, the metric we use 
is a reasonable approximation, and the resulting entanglement entropy should be similar. 
In bulk language, we expect that as long as other fields are not turned on too strong and so 
do not backreact too much on the metric, this result should be a very good approximation. 

Note that in the region Ilia the entanglement entropy agrees with the entanglement 
entropy for an extremal BTZ black hole. One can extend our work to easily show that if we 



have left- and right-moving packets and we measure the entanglement entropy using (3.7) 



in the overlap region we find the entanglement agrees with the nonextreme BTZ black hole 
0, 



S. 



overlap 



6 



log 



2 sinh[r + (i] sinh[r_rf] 



T. 



±± 



8ttG 



N 



(3.13) 



T + T_<5 2 

However, the overlap region does not extend to late times, which means that the system 
does not actaully thermalize. This is because we are only considering vacuum gravity in 
2+1, and so the left and right moving sectors do not interact. 
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A=2,t=1/8, d=l, 3 



A=2, t=10, d=l, 3 





x+ 



FIG. 5: Plots of AS = S — | log [2<i 2 /5 2 ] at fixed r and A as we increase d. The lightest curve 
has d = 1 < A and the darkest curve has d = 3 > A. 



Thanks to the fact that we know the geometry is exact AdS, we can use the arguments of 
[5] to know that we are in fact calculating the covariantly defined holographic entanglement 
entropy It is worth pointing out that the only other possible geodesic connecting the two 
points is a piecewise smooth geodesic which falls down to the horizon from the two points, 



and then runs along it. We know from the form of our metric (2.1) that the horizon sits 
at z 4 = 4(L + L_) 1 / 4 , and, thanks to the fact that we only have left-moving stress tensor 
turned on the horizon sits at z = and a geodesic falling down will always have an infrared 
divergent length, and therefore the shortest geodesic is the connected one whose length we 
have used. 



IV. SHOCKWAVE LIMIT 

Now we will discuss what happens when we try to make our pulse of energy look more 
like the traditional AdS Aichelburg-Sexl shockwave |1T] (for a review see eg [13]). Recall 
that the shockwave metric in AdSd+i is 

i 2 

ds 2 = — \-2dx+d- + dx% + dz 2 ] + $(x T , z)5(x + )dx 2 + (4.1) 

z 

where $ is a nontrivial function. We see that our piecewise geometry is a simple example 
of a shockwave, but since we do not have transverse directions in AdS% we do not have 
to support all of the stress tensor on a delta function. Noting that the total energy on a 
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constant t slice is 
E 



T 00 dx 



167rG 



N 



t + X 

2V2 



dx 



16nG 



-V2Ar 2 



(4.2) 



A? 



Keeping e = At 2 finite while sending A — > simplifies things, as now the pulse is localized 



solely at x + = 0. Comparing (4.1) to our metric we find $ = £ 2 \/2e, which agrees with the 



d = 2 shockwave metric. The entanglement entropy reduces to 



AS 







log 



1 + 



EX+(d— £+) 







x + < 
< x + , x + — d < 



(4.3) 



Note that in this limit the complicated behavior in the entanglement has simplified dramat- 
ically, and we only have an entanglement entropy different from the vacuum value when the 



pulse is within the region A. The maximum entanglement entropy (3.12), now occuring at 
= d/2, has simplified to 

maxfAS] = % log[l + del A], (4.4) 
6 

which features logarithmic scaling for large de. Since L + is now a delta-function and not 
a piecewise-continuous, the entanglement entropy has a discontinuous first derivative. In 
figure [6] we provide plots of AS as we send A — > 0, with fixed e. 

d=2, e=l, A=2,...,0 d=2, e=100, A=2,...,0 



0.05 
AS °.04 
c 003 
0.02 
0.01 
0.1)0 





1.5 

X+ 



1.5 



FIG. 6: Plots of AS at fixed d and e as we send A 
lightens as we approach A = 0. 



0. The darkest curve is A = 2 and it 



This result is very similar to CFT calculations of entanglement entropy of 1+1 systems 
subjected to local quenches, which are constructed by adding or removing a point-like defect 
seperating a spacial slice into two regions [6J [7] . Changing the topology of the spacetime 
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will cause a local UV divergence in the form of an infinite-energy pulse leaving the defect at 
the speed of light. If we join two manifolds separated by a co-dimension one boundary with 
transverse translational symmetry, it is reasonable to expect that this will source a bulk 
gravitational shockwave, which we can model as an Aichelburg-Sexl shockwave. Indeed, it 



the limit e ^> 1, we find (4.3) looks very similar to the CFT results, where we identify e with 
the UV regulator of the local quench, and we find that there is no residual entanglement 
entropy after the pulse has passed, indicating we have lost information of the boundary 
entanglement entropy. This one missing component can be included by assuming there is a 
contribution to the initial entanglement entropy coming from the co-dimension one boundary 
defect, treating it perhaps along the lines of [14"] . 

The entanglement entropy for the AdS% shockwave can be generalized to higher dimen- 
sions. Further we can use the fact that it can be constructed by taking a singular limit of 
a boosted global AdS black hole. We can then use the covariant holographic entanglement 
entropy construction and apply it to nontrivial time-depentent boundary regions and find 



a higher- dimensional analog of (4.3). While we leave this calculation to future work, we 
conjecture that in general the form will be similar, with maxfAS 1 ] oc (de) D ~ 2 as de 3> 1. 
where we are working with a D-dimensional CFT, d is the length scale of the boundary 
region A, and e is the energy of the shockwave. 



V. CFT COMPARISON 



Calculating an entanglement or Renyi entropy in a generic state in a CFT is quite non- 
trivial. However when we consider states constructed by acting with a nontrivial conformal 
transformation on the vacuum, we can use the fact that the replica trick involves twist op- 
erators which transform as primaries under conformal transformations |15j . This argument 
has been used previously to calculate the entanglement entropy for a thermal state or the 
theory on a circle. If we can analytically continue the argument of [15] to Lorentzian signa- 
ture, we expect that trp^ transforms under scale and conformal transformations as the nth 
power of a primary operator O n with A + = A_ = £(1 — 1/n 2 ) = A n , where A± are now 
left and right central charges. This is a nontrivial assumption, because the formula for trp^ 
is derived from a Euclidean path integral, and we are considering states which can not be 
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represented as Euclidean path integrals. In vacuum, for a single interval, 



trp* 



2(y 2 + -yi+){yi- 



V2- 



where c n are not generically determined, except for c\ = 
lattice scale. Recall that under the map x± y± — f±(x± 



= (O n ( yi )O n (y 2 )) n . (5.1) 
1, and a is a nonuniversal UV 



{O n ( Xl )O n (x 2 )) n = |/i(x 1+ )/i(x 1 _)/i(x 2+ )/i(x 2 _) 



inA; 



2 If + (X 2+ ) - U(X 1+ )} - /_(X2_)] 



(5.2) 



We can now evaluate Sa 



S 



d n [trp^] | n=1 and we find exactly (3.7), 
2A/ + (x 1)2 )A/_(x lj2 ) 



6 



log 



'V /+2/+1/-1/- 



(5.3) 



We wish to emphasize that since this is so far just an analytic continuation of the behavior 



of euclidean path integrals, (5.3) is only a conjecture for generic coordinate transformations. 
While one would like to simply say this is a consequence of conformal invariance of the 
theory, the statement is not quite so straight-forward. The important f' ± contribution to 
the entanglement came from very different places on either side of the duality. On the field 
theory side, it is a consequence of the fact that the twist operators transformed as primaries. 
However, on the gravitational side, it came about by carefully treating the transformation of 
the UV regulation of the geodesic length. The simplicity with which this worked out suggests 
that transformation laws for more complicated objects, such as mutual Renyi information, 
or mutual entanglement information, for multiple disconnected regions may have similar 
transformation properties. It would be very interesting to better understand holographic 
Renyi entropies to verify that they satisfy the transformation properties given above. 



VI. CONCLUSIONS 



In this work, we have calculated the time evolution of the entanglement entropy of a 1+1 
CFT subjected to a localized packet of energy density. As a tool in doing this calculation 
we have integrated the infinitesimal coordinate transformations that shift the stress tensor 
of the CFT, providing a closed-form diffeomorphism that takes AdS^ in the Poincare patch 
to the most general solution. We have also verified that our result agrees with a CFT 
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calculation of entanglement entropy of a single interval for a state created by acting on the 
vacuum with a coordinate transformation. We find that in the narrow shockwave limit, 
the form of the entanglement entropy simplifies greatly, and reproduces some but not all 
of the structure of CFT calculations of time evolution of entanglement entropy using CFT 
techniques [HI [7] . 

There are many interesting future directions one may wish to look towards related to 
this work. First of all, the ability to construct nontrivial nonperturbative time- dependent 
backgrounds in asymptotically AdS 3 is a powerful tool. It would be interesting to study 



general retarded correlators on (2.1) by mapping general solutions to the Poincare patch 
slicing. We can use the WKB approximation to relate the bulk geodesic length between two 
space-like separated points to a two-point function for high dimension operators [EJH7], so 
identifying G(x\\X2) ~ e~ mxLength ( 7a: i' :E 2) with m = A/£ + . . ., we can generalize (3.11). In 
the shockwave limit we find 



(2Ax+Ax_) A | x+ + Ax + < 



G(x±\x±T Aar±) = <T A < (2(Ax+ - x + (x + + Ax+)e:)Ax_) | x + < < x + + Ax+ 

(2Ax+Ax_) _A | < x + 

(6.1) 

where we have assumed that Ax± > 0, ensuring that the points are spacelike separated. 
It would be very interesting to extend this result to the two-point function for a classical 
field in the shockwave spacetime. Similarly it would be very interesting to use the singular 



diffeomorphism which gives us the higher dimensional shockwave metric (4.1) and study 
nontrivial time-dependence of correlators there. 

It would also be very interesting to study various methods of quenching a holographic 
system which would source the Shockwaves we consider, such as inserting or removing a 
holographic defect [TH [19] or a boundary [TJ]. It would also be interesting to study the 
inclusion of nontrivial behavior of other primary operators either in the field theory directly 
or as additional bulk fields, and see if with their inclusion the system could thermalize. 
Whether we could continue to use purely analytic techniques or would need to resort to 
numerical methods is unclear. 
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